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invariant mass terms. 
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In very special relativity (VSR) Q], the laws of nature are not invariant under the full 
Lorentz group but only a subgroup called SIM(2). In SIM(2) invariant theories there is 
a preferred light-like four- vector, called n fl . The unbroken generators of SIM (2) consist 
of the generators of little group of n^ 1 as well as boosts under which transforms as 

— > e a n^. SIM(2) invariance is sufficient to guarantee that particles obey the relativistic 
dispersion relation, E 2 = p 2 + m 2 . This ensures that SIM(2) invariant theories will satisfy 
constraints that come from relativistic kinematics, so the violation of Lorentz invariance in 
these theories is more subtle. Since SIM(2) has only one-dimensional representations, states 
which are degenerate in a Lorentz invariant theory can split into nondegenerate states which 
satisfy relativistic dispersion relations with different masses. For instance, we will see below 
that for certain kinds of SIM(2) invariant mass terms, the four spin states of the electron 
and positron split into two nondegenerate pairs of states. SIM(2) theories respect CPT [lj, 
which accounts for the remaining degeneracy. The other way the SIM(2) invariant theories 
can violate Lorentz invariance is via interactions. For example, the couplings of an electron 
to background electromagnetic fields, particle decay rates, and cross sections can depend on 
the preferred four-vector and therefore violate Lorentz symmetry. These effects cannot 
be accessed simply by examining free particle propagation. 

When the SIM(2) generators are combined with parity operation, the full Lorentz group 
is obtained. Therefore, the breaking of Lorentz invariance is tied to the breaking of parity. 
This may naturally suppress Lorentz violating effects. Another important feature of the 
SIM(2) group „ that it has no invariant tensors, implying that the nsna. spurion analysis 
of Lorentz violation does not apply to SIM(2) invariant theories 111. In traditional analyses 

nnnn 

of Lorentz violation, see e.g., Refs. [2J, |3|, |5( , Lorentz noninvariant tensors are introduced 
and then local dimension-3 and dimension-4 operators in the Standard Model containing 
these tensors are added to the Standard Model Lagrangian. This approach does not lead to 
SIM(2) invariant theories, since the terms with Lorentz noninvariant tensors will also violate 
SIM(2) invariance. In SIM(2) invariant theories, Lorentz violating terms in the Lagrangian 
are necessarily nonlocal. 

Finally, SIM (2) invariant theories provide a novel way of generating a mass for a two- 
component spinor (f| . This is especially interesting for theories of the neutrino mass. To give 
neutrinos a mass in the Standard Model one must either add a gauge singlet right-handed 
neutrino to form a Dirac mass term or else add a Majorana mass term which violates lepton 
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number. SIM(2) invariant theories allow for a novel mass term which requires neither a 
right-handed neutrino nor lepton number violating interactions. Ref. p proposed a SIM (2) 
invariant mass term for neutrinos and studied its effect on the electron spectrum in tritium 
beta-decay. 

When the SIM(2) invariant mass term of Ref. p is made invariant under the SU(2) l x 
U(l) gauge group of the Standard Model, Lorentz violating interactions will be induced in 
other sectors of the theory, including QED. The purpose of this paper is to study the Lorentz 
violating effects implied by the SIM(2) invariant neutrino mass. We also make a preliminary 
investigation of loop effects in a SIM(2) invariant extension of QED. 

The paper is organized as follows. In the section I, we generalize the mass term for neu- 
trinos considered in Ref. p by requiring the theory to be invariant under the SU (2)^ x U (1) 
gauge symmetry of the Standard Model. This leads to Lorentz violating effects in QED as 
well as tree-level lepton family number violating interactions. In section II, we focus on a C 
invariant and SIM(2) invariant version of QED, called VSRQED, and derive a nonrelativistic 
Hamiltonian for electrons in the presence of background electric and magnetic fields. This 
Hamiltonian is used to determine the consequences of VSRQED for tests of Lorentz invari- 
ance with trapped electrons. VSRQED corrections to the spectroscopy of hydrogenic atoms 
are also considered. Since Lorentz violating effects in this theory are suppressed by m1/m 2 el 
the predicted Lorentz violating effects are tiny and within current experimental bounds. In 
section III, we study one-loop self-energy corrections to the photon and electron propaga- 
tors. The nonlocality of the SIM (2) invariant theory leads to IR divergences which make the 
loop integrals ill-defined. We propose a prescription for dealing with these IR divergences. 
Evaluated with this prescription, the loop diagrams give corrections that satisfy the Ward 
Identity. In addition, the photon remains massless at one loop and no Lorentz violating 
modifications to the photon propagator are induced. Furthermore, the one-loop correction 
to the SIM(2) invariant electron mass term is finite. In section IV, we calculate the tree 
level contribution to the lepton family number violating decay fi —>■ e + 7. Our results are 
summarized in section V. In Appendix A, we derive Feynman rules for VSRQED. 
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I. SIM(2) INVARIANT MODIFICATIONS OF THE STANDARD MODEL 



The SIM (2) invariant, Lorentz violating, modified Dirac equation for the left-handed 



neutrino proposed in Ref. 



.1 



is 



j j-^JL)^(p) = o. (i) 



2 n ■ p / 

Though Eq. violates Lorentz invariance, by squaring the modified Dirac operator it is 
obvious that 

(p 2 -ml)v L (p) = 0, (2) 

so the left-handed neutrino has a relativistic, massive dispersion relation. Thus, SIM(2) 
invariance provides an alternative way to introduce a neutrino mass without adding either 
a Dirac or Majorana mass term. 

The Lagrangian that would give rise to Eq. is 

£ = ^-™>- < 3 > 

where vl is the left-handed neutrino field. The second term in the Lagrangian is nonlocal 
and has a structure similar to kinetic terms in th e lig ht-cone QCD Hamiltonian [?} Is| as 
well as Soft-Collinear Effective Theory (SCET) P, Il0| . The action is ill-defined for modes 
of vl satisfying in ■ Ovl = 0. In the tree level calculations of this paper we will not have to 
worry about this issue, however, in loop interactions this leads to infrared divergences and 
a prescription is needed to define the theory. We postpone discussion of this issue until we 
study loop diagrams in section IIHI 

In this section we focus on the implications of embedding the Lagrangian of Eq. (jSJ) 
in the Standard Model. A necessary requirement is that the action be invariant under the 
SU (2)lxU(1) gauge group. Eq. (jH)) can be made gauge invariant by replacing the neutrino in 
Eq. Q with the left-handed lepton doublet and making the derivatives in Eq. © covariant. 
We can also consider a similar term for the right-handed electron field. We will restrict our 
attention to the lepton sector of the Standard Model. For three generations of leptons, we 
obtain the following SIM(2) invariant mass terms, 

= ~T» Mi - \e R „ Ml jL^ , (4) 
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where the SU(2) L doublets are denoted by II, % and j are family indices, and M 2 and M 2 
are Hermitian matrices. Summation over repeated indices is implied. Once the electroweak 
symmetry is broken the mass terms in the theory are 

1 ji 
C = -v(Xij e Lti e RJ + h.c.) - -P £ji M 2 - —j^Lj 

where ei,i and e R ^ are the left-handed and right-handed charged lepton fields, respectively, v 
is the Higgs expectation value, and the are the usual Standard Model Yukawa couplings. 
Three independent unitary transformations in family space, 

eL^U L e L , e R ~^U R e R , v L -> V v L , (6) 

can be used to diagonalize some of the mass terms. The matrix Vm = U\V is the neutrino 
mixing matrix. Not all SIM (2) invariant mass terms can be simultaneously diagonalized. 
Choosing U L and U R to diagonalize the charged lepton Dirac mass matrix and V to diago- 
nalize the neutrino mass matrix, we find 

1 ,-. ™2 i 



£ = -m e ,i (e Lji e R>i + h.c.) - -u Lji m uii -—^vl 



~^e Ltl (V M ml V^ij — ^ e LJ - ^e R>l (U R M 2 U R ) l3 — ^ e RJ , (7) 

where m^i would be the charged lepton masses in the absence of SIM(2) violating terms, 
and we have defined the diagonal matrix of neutrino masses squared to be m 2 u = V^M 2 V. 
The SIM (2) invariant mass matrix for the left-handed leptons is determined in terms of the 
neutrino masses and mixing matrix. It is nondiagonal and therefore these terms introduce 
additional lepton family number violation. The SIM(2) invariant mass matrix for the right- 
handed leptons is undetermined and should generically be nondiagonal since the matrix U R 
which is needed to diagonalize the Yukawa couplings, A^, is in general not the same as the 
one required to diagonalize Mfj. 

The existence of different SIM (2) invariant mass matrices for the left-handed and right- 
handed leptons leads to a wave equation for charged leptons that is more complicated than 
the neutrino's wave equation. In particular it can lead to a breaking of the degeneracy of 
the spin states of the charged leptons. We will focus on the electron since we are primarily 
interested in looking for Lorentz violating effects in SIM (2) invariant extensions of QED. 



Lorentz violating effects which split the degeneracy of the electron's spin states along a 
preferred direction in space are strongly constrained by modern experiments with trapped 
electrons 

For the electron the modified Dirac equation takes the form 

^ - m - \^i M l + M -l^j «(p) = , (8) 

where u(p) is electron spinor, m = m^i, and M± = ±-(^777^^)11 + ^(U^^UrJu. Note 
that the combination M 2 — M 2 = (yjfm^V^)n is the so-called electron based neutrino 
mass squared, mi. The orthogonal linear combination M 2 + M 2 = (Up l M 2 Uii)ii is not 
determined by parameters in the neutrino sector. 

The current constraint on mp from tritium beta decay experiments is rap < 2eV jl^ . 
This constraint assumes a Lorentz invariant neutrino mass. Ref. 0] pointed out that the 
effect of a SIM (2) invariant neutrino mass on the beta decay spectrum is roughly equivalent to 
that of a Lorentz invariant neutrino mass that is a factor of two more massive. So for SIM (2) 
invariant masses, the current constraint from tritium beta decay should be interpreted as 
mj3 < 1 eV. Cosmology gives a comparable constraint on the sum of neutrino masses, which 
must be greater than mp. The bounds are somewhat model dependent, with quoted lower 
bounds ranging from 0.6- 2.0 eV Q,Q,Q- 

We do not know how the cosmology bounds 
would be affected if the neutrino masses are SIM(2) invariant rather than Lorentz invariant. 
We will assume mp < 1 eV for the remainder of the paper. 

To determine the mass spectrum, we invert the operator in Eq. (JHJ) and find the locations 

of the poles in the electron propagator. First we define 

M 2 

P>^P»-^—, (9) 
2 n ■ p 

and note that n-p = n-p' and p' 2 = p 2 — M 2 ^. Left multiplication of the operator in Eq. (JHJ) 
by f - m + %^75 yields 

p' 2 + m 2 - 2m p" - M! 75 . (10) 
Next, multiplying by 2m ^' + {p' 2 + ttiq) + M 2 ^ 5 yields the operator 

(p'2 _ m 2 ) 2 - Mt . (11) 
This sequence of operations shows that the electron propagator is 

, (2m4 - m M 2 + ^- p + (p 2 + m 2 - M 2 ) + M 2 ^ - m - \j- p (M 2 + - M 2 _ l5 )) 

% (p2 _ m 2 _ M 2 y _ M 4 > ( 12 ) 
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so the electron propagator will have poles at 



V 



2 - E 2 ± = ml + Ml ± Ml 



(13) 



The four spin states of the electron and positron are degenerate only when M_ = 0. 

To better understand how the states are split we study the solutions to the modified 
Dirac equation. We work in the rest frame p M = (E,0), and choose the z-axis to be the 
spatial direction singled out by the four-vector, n M = (1,0,0,1). Then the modified Dirac 
equation is 



{ El °- m °-^E {l ° 



7 3 )(M2+M! 75 )U(p)=0 



(14) 



Using the Bjorken-Drell conventions for the gamma matrices, we find the following solutions: 

/ \ 
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where E± = -J m\ + M\ ± M£ . We have labeled the positive energy spinors u-j- and ttj_, and 
the negative energy spinors v\ and Vy. In the absence of SIM(2) invariant mass terms, u\ (itf ) 
corresponds to an electron with spin-down (spin-up) along the z-axis. The negative energy 
solutions are interpreted as antiparticle states which are degenerate with the particle states, 
as expected because the SIM(2) invariant masses respect CPT. We see that when M_ ^ 
electrons with spin- up and spin-down along the z- axis are split by E + — E_ w Ml/(2mo). 

Such a splitting will lead to a diurnal variation of the anomalous precession frequency of 
capped electrons in expe„ m en ts that m e„ - 2. Ref. Q place, bonnds on dinrnal 
variation of the anomalous precession. Ref. [HJ analyzes the data using the formalism of 
Ref. 16], which considers an electron wave equation of the form: 



{j)-m e - eA^ - - b^ 5 ^)u(p) = 



(15) 



7 
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P 


T 






#i^^(*,-x) 


#^a^(-i,x) 


#75 1^3 V^*) 


-#75^^,x) 


-#75^g^(* 5 -x) 





TABLE I: Action of discrete symmetries on SIM(2) invariant mass terms. 

where is the electromagnetic vector potential, and and are background four- vectors. 
The vector has no physical effects since it can be removed by the gauge transformation 
— > A^ — d^A, where A = a^x^/e. The four- vector 6 M gives rise to a diurnal variation of the 



anomalous precession frequency (ill . . Eq. ()15j) does not have the same form as Eq. (JBJ) . 
Invariance under n M — > e a n' 1 allows us to rescale so that n 11 = (1, n), where n is a spatial 
unit vector. Then for nonrelativistic electrons we can approximate n ■ y « m e and identify 



Mi/(2m e ) n„ and 6^ = M 2 /(2m e ) n„. The observation of Ref. 



that does not 



affect the precession frequency is consistent with our observation that the term proportional 
to M\ does not split the electro n sp in states, and the diurnal variation of the anomalous 
precession due to found in Ref. [la] is consistent with the splitting of electron spin states we 
find in the presence of the term proportional to M 2 . Ref. obtains a bound on a product 
of a certain component of 6^ times geometrical factors associated with the orientation of 
the experiment relative to the background vector. Assuming these geometrical factors are 
of order unity, we find M 2 _ = -l(V M mlV^)n + l(U j R M 2 U R )n < l(r 10 eV 2 . This value of 
M_ appears to require considerable fine tuning since (Vm^V^-)ii = mi > 1CT 5 eV 2 , as we 
will later show. 

However, it is possible to enforce M_ = by imposing charge conjugation invariance 
on the SIM(2) invariant theory. The action of discrete symmetries on the terms in the 
Lagrangian which give rise to the SIM(2) invariant mass terms of Eq. (|SJ) is given in Table UJ 
Here n M is the parity image of n M , i.e., if n M = (1, n) then n M = (1, — n). From TableQwe see 
that both SIM(2) invariant mass terms respect CPT, which is consistent with the arguments 
of Ref. The mass term without 75 is also invariant under C and PT, while the mass term 
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with 75 is not invariant under either. Thus, if we impose C invariance, or equivalently PT 
invariance, we can get rid of terms proportional to M£. For the case = 0, the electron 
propagator is 




which is considerably simpler than the electron propagator in Eq. (|12|). Thus, imposing C 
invariance makes the SIM (2) invariant version of QED easier to analyze. 

There are two important conclusions we wish to draw from the analysis of this section. 
One is that SIM(2) invariant neutrino masses necessarily require SIM(2) invariant mass terms 
for the electron, implying Lorentz violating effects in QED. Tests of Lorentz invariance in 
QED have reached a high level of precision so the hypothesis that the neutrino has a SIM (2) 
invariant mass term can be tested by looking for Lorentz violating effects in this sector of 
the Standard Model. Current experiments with trapped electrons force one to consider a C 
invariant and SIM(2) invariant extension of QED, which we will call VSRQED. This theory is 
studied in the next two sections of the paper. The second important point is that the SIM (2) 
invariant mass terms in the lepton sector are family non-diagonal so SIM(2) invariant mass 
terms for neutrinos imply new tree-level contributions to lepton family number violating 
processes such as \i — > e + 7. This process is studied in section IV. 

II. NON-RELATIVISTIC ELECTRONS IN VSRQED 

In this section we derive the coupling of nonrelativistic electrons to electric and magnetic 
fields and study the effect these interactions have on measurements of g e — 2 using trapped 
atoms and on the spectroscopy of hydrogenic atoms. We will impose C invariance so the 
modified Dirac equation for the electron is 



where u n (p) is the electron spinor and m/| = (VmitiIV^ii is the electron based neutrino 
mass. For the remainder of this paper, u n (p) will denote a positive energy solution to the 
modified Dirac equation in Eq. (|17j) . and u(p) will denote will denote a positive energy 
solution to the unmodified Dirac equation. It is useful to write the solution u n (p) in terms 




(17) 
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of u(p), 

f m e -m i \ 
u n (p) = [1 - — ) u(p) 

1-^^W (18) 
4m e n ■ p J 

Note that m e = yfm% + m|, — m e )u(p) = 0, and in the last line we have expanded to 
first order in m 2 ^/m 2 . This allows us to write the interaction between the electron and the 
electromagnetic fields to 0{rni/rn 2 e ) in terms of solutions to the conventional Dirac equation. 
If we work in n ■ A = gauge, the Feynman rule for the interaction between the electron 
and the photon is 

-ieu n (p')ru n (p)A,(q) = -ieu(p')ru(p)A,(q) + i--J-u{p')^- f f. u ( p ) + ... . 

4 77i e \n • p ) \ n ■ pj 

For nonrelativistic electrons, we can expand to lowest order in the electron momentum and 
express the interaction in terms of two-component non-relativistic spinors. To lowest order 
in the electron momentum, we find 

.em} n v n a a aii F^{q) em} 

*t— u[p)— ^ T-u[p) « ^ V T (( n x E) • a - a • B + (n • cr)(n • B))ip . 

4 m e [n ■ p ) [n ■ p) 2 mr e 

We have set = (1, n), where n is a unit vector. This leads to the following contributions 
to the Hamiltonian for a nonrelativistic electron in the presence of E and B fields: 

H S im{2) = e/2 B [-(n x E) • a + o • B - (n ■ <r)(n • B)] , (19) 

where \xb = eh/(2m e c) and e = m^/m 2 . Currently, e is constrained to lie within the range: 

9.2 1(T 17 < e < 3.8 1(T 12 . (20) 

The upper bound comes from mp < leV [12], and the lower bound comes from noting that 

2 /t t 2-rzt \ 22 2,22 2,2 2 

m p = [V M m v Vh)u = c 12 c lz m t + s 12 c 13 m 2 + s 13 m 3 

= m\ + s\ 2 c 2 l3 /\m 2 2l + s 2 13 Am 2 31 . (21) 



Here c^- = cos(0y), s^- = sin(6 l i j) and Am| = m 2 — m 2 . This is a sum of positive definite 

0. 

2.410" 5 eV z . There are no lower bounds on m{ or s 13 . 



terms. With Am^ = 8.0 10~ 5 eV 2 , s\ 2 = 0.31 and c 2 13 > 0.95, we find mj > s 2 2 c 2 3 Am^ 
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For nonrelativistic electrons the SIM (2) invariant mass term gives rise to a Lorentz invari- 
ant contribution to g e — 2 as well as Lorentz violating couplings of the electron to background 
electric and magnetic fields. These interactions depend on the preferred direction, n. Note 
that manifest invariance under the SIM(2) transformation — > e a n^ is lost when the non- 
relativistic limit is taken. The unit vector, n, points in the preferred direction as seen from 
the electron rest frame. The size of the Lorentz violating corrections are fixed by parame- 
ters in the neutrino sector. Therefore, the hypothesis that neutrino mass terms are SIM(2) 
invariant rather then Lorentz invariant predicts tiny Lorentz violating effects in QED whose 
size is determined by e. Tests of Lorentz and rotational invariance in QED which constrain 
e to be less than 9 10~ 17 will rule out the SIM(2) invariant neutrino mass hypothesis. Ob- 
servation of Lorentz violation in the range predicted by Eq. (119)) and Eq. (|2(J|) would provide 
support for SIM(2) invariant neutrino masses. We will consider the effect of the terms in 
Eq. (|19j) on experiments that search for an electric dipole moment (edm) of the electron, 
experiments with trapped electrons, and the spectroscopy of hydrogenic atoms. 

The first Lorentz violating term in Eq. (fTH|) is like an electron edm, except the electron 
spin couples to n x E instead of E. Searches for an electron edm use an apparatus in which 
atomic beams are passed through regions in which E and B fields are either aligned or 
antialigned, and then an edm will result in a splitting between the atomic levels in the two 
beams, which can be measured via interference 17]. These experiments are insensitive to 



the n x E • a interaction in Eq. ()19)) because the interaction vanishes when E and B are 
aligned. We are not aware of experiments which look for splittings between spin states of 
atoms in crossed E and B fields, which is what would be required to test for interactions 
proportional to nx E-cr. Note that the current bound on the electron edm is \d e \ < 1.6 10 -27 
e-cm = 8.3 10~ 17 \xb- Thus a bound on e/ig comparable in size to present bounds on the 
electron edm requires e smaller than the lower bound in Eq. (|20|) and would rule out the 
SIM(2) invariant neutrino mass hypothesis. 

The second term in Eq. (fTH|) is a non-Lorentz violating correction to the electron 
magnetic dipole moment. The most recent experimental measurement in Ref. ^| finds 
g e /2 = 1.00115965218085(76). This is consistent with the theoretical calculation of Ref. ^| 
which finds g e /2 = 1.00115965217586 ± 8.48 10~ 12 . The second term in Eq. (0 gives a cor- 
rection S(g e /2) = e, which leads to the bound e < 1.27 10~ n . This is slightly larger than 
the upper limit in Eq. (|20|). 
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The third term in Eq. (fT^j) modifies the precession of the electron's spin in a background 
magnetic field. The effect of this term is to shift the effective magnetic field seen by the 
electron spin: 

H = -y/UB a B' = -^Hb o- ■ B - e/i B (n • a) (n • B) , 

where B' = B + (2/g e ) en (n • B), so the the modified precession frequency is 

u = g e fi B \B , \ = g e fi B \B\ (\ + |-e(n-B) 2 ^ . (22) 

This will result in a time dependent precession frequency since B is a magnetic field in a 
rotating frame on the surface of the Earth, while n is a constant vector which is assumed 
to be fixed with respect to distant stars. Assuming the B field is oriented perpendicular to 
the surface of the Earth and outward, the magnitude of the daily variation in the precession 
frequency is 



UJ 



2 I sin(2A) sin(2a)| |tanAcota|>l 

— e < V ' K n ' ~ . (23) 

9e 1 max[sm 2 (X ± a)] | tan A cot a\ < 1 



Here u is the time-independent precession frequency in the absence of the third term in 
Eq. (|19jh Au; is the maximal daily variation in u, A is the latitude of the laboratory exper- 
iment, and a defines the angle between n and the axis of rotation of the Earth. The bound 



on Aw from Ref. [ll| is Aw < 1.2 Hz. The experiment is conducted with a B-field of 5.85 
T, corresponding to a precession frequency of ujq = 5.1 10 11 Hz. Assuming the geometrical 
factors in Eq. (J23)) are order unity we find e < 10 -11 , which also exceeds the upper limit in 
Eq. (|2Uj). 

Next, we consider the effect of these terms on the atomic spectra of hydrogenic atoms 
applying perturbation theory in e. The factor of (n x E) • S becomes Ze n • (r x S)/r 3 when 
the electric field of the atomic nucleus is substituted in for E. In this form, we see that this 
term does not shift the energy spectrum to first order in e as the matrix element of this 
operator vanishes between any two states with the same n. We have 

r 



[p, H ] = -iZe" 
where H Q is the Coulomb Hamiltonian, so 



^e n " ( 7 S) = V(pxS),g,]. (24) 



r 3 e ' 
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Therefore 



(n,l\Ze n '^ S \ n, I') = {n,l\ % -[xi- (p x S),H Q ]\n,l') = 0. (25) 

The third new term in the non-relativistic Hamiltonian gives a rotationally noninvariant 
spin-orbit coupling because of the magnetic field in the electron's rest frame generated by 
the rotating nuclear electric field. The correction to the Hamiltonian is 

A# = e-4-L-nS-n 

2wig r 6 

The directional nature of this term causes a splitting in originally degenerate 2p 3 / 2 states. 

Choosing the n direction as the axis of quantization, and using 

/ , i 1 i , \ Z 3 a 3 ml 

(n, I, m h m s \ — \n, I, m u m s ) - 



we see that the originally degenerate energy levels split by 

Z 4 a*m e ZW(Ry) 

Hi r l 3 — rLi r i i = 6 — 6 

J z -±- J z _±- 72 36 

This splitting must be less than the precision of the n_nt of the 2s 1/2 - 2p 3/2 interval 
in hydrogen measured in Ref. 20], otherwise two distinct lines would have been observed. 
Using the measured splitting frequency 9911.200 ± .012 MHz, we find e < 2.5 • 10~ 6 . 

Furthermore, this term causes a correction to the n = 2 Lamb shift: it leaves the 2si/2 
states unshifted while shifting the 2pi/ 2 states. Using 

2p,j = ^rrij = ^\ = \mi = 0,+) - ^| \mi = +1,-) 

2p,j = ~,mj = ~\ = = -!>+) - \j\\ m i = °'~) 

We find that the energy of 2pi/ 2 states shift by 

ZW(Ry) 
Ai/ 2pi/2 - -e — 

This is a correction to the measured Lamb shift, |Aa;£ 0TO &| = 1057.839 ± .012 MHz j^], 
which sets an upper bound for e at 4.9 10~ 6 . 

To summarize, bounds on Lorentz violating and Lorentz preserving couplings to magnetic 
fields from measurements of g e /2 with trapped electrons yield e < 8.5 10 -12 , which is just 
above the upper limit in Eq. fl2()H . Bounds from the spectroscopy of hydrogenic atoms are 
significantly weaker, with e < 10~ 6 . Bounds on a Lorentz violating electric dipole moment 
comparable to existing bounds on electron edm can rule out SIM(2) invariant masses. 
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III. ONE LOOP SELF-ENERGIES 



In this section we consider the one-loop corrections to the photon and electron propaga- 
tors. Many high precision tests of Lorentz invariance involve the propagation of light, for 
example, tests of cosmic birefringence. Therefore we would like to know whether Lorentz 
violating modifications of the photon propagator arise in SIM(2) invariant theories. Since 
we have not added a SIM (2) invariant modification to the gauge field kinetic terms, such 
corrections must occur via loop effects in this theory. 1 It is also important to understand 
how the nonlocality of the SIM (2) invariant action affects the calculation of loops. Straight- 
forward evaluation of one loop corrections to the photon and electron propagators results in 
unregulated infrared (IR) divergences. These are due to the factor of l/(n -p) in the electron 
propagator. These IR divergences cannot be regulated by dimensional regularization or by 
taking external legs in the diagrams off-shell, as is the case for IR divergences that arise 
in local quantum field theory diagrams. This is not so surprising given the nonlocality of 
the fermion action. It indicates that the theory is ill-defined at the quantum level, unless a 
prescription for dealing with these IR singularities is given. 

Since the physical origin of the nonlocality is unknown, it is not clear how to construct an 
adequate prescription. Motivated by the observation that the action is ill-defined for off-shell 
modes with n ■ p = 0, we suggest that these modes be excluded from the theory altogether. 
Since on-shell momenta cannot satisfy n ■ p = 0, these modes only appear in internal lines 
in Feynman diagrams. We propose that in loop diagrams n ■ p = modes should not be 
allowed to propagate. This can be implemented by imposing a "+" prescription for defining 
the quantum mechanical theory, in which loop integrals over the n ■ I component of loop 



1 It is possible to add a SIM(2) invariant modification of the gauge field kinetic term |2l| : 

I—L — _ — n a F afl — - — TigF^ . (26) 
2 in-D in-D P » y ' 

In n ■ A = gau ge this is simply a photon mass term. The present bound on the photon mass is 
m 7 < 610 -17 eV |l2j . which is about 15 orders of magnitude smaller than mass scale set by neutrino 
mass differences. For this reason we have chosen to ignore such a term. 
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momenta I are defined to be 



n ■ I J [n ■ l) + 

tn.l'W-M (27) 

This prescription tames the IR singularities. The novel contribution to the photon self 
energy in the SIM (2) invariant theory vanishes on-shell with this prescription. At one-loop, 
the photon propagator is essentially unmodified relative to the usual self-energy correction 
in QED. There are no Lorentz violating effects and no photon mass is induced. This last fact 
is important since it is possible to write down a SIM(2) invariant, gauge invariant photon 
mass term (see footnote). If an 0{ocmp) mass term were induced by radiative corrections 
the theory would be phenomenologically unviable. The "+" prescription also removes IR 
divergences in the electron self-energy. Moreover, in that case we find that the radiative 
correction to the SIM (2) invariant mass term is finite. This is something that one might 
expect from the naive degree of divergence of the SIM(2) invariant correction to the electron 
propagator. Counting the naive degree of divergence in higher order diagrams suggests that 
the finiteness of the correction to the SIM(2) invariant mass term will persist in higher orders 
in perturbation theory. If this is the case, then SIM(2) masses are natural in the technical 
sense. To have small SIM(2) invariant mass terms, presumably of order the neutrino mass, 
will not require fine tuning at each order in perturbation theory. However, the question of 
why the scale of SIM(2) invariant masses is so much smaller than the other mass scales in 
nature remains unanswered. 

We begin by computing the electron self-energy in n ■ A = gauge. The diagram is the 
same as in QED and the result for one-loop self-energy is 

( _ , m % j\ 

_TM 2 f d D g ^\ i + m »+ 2 n- q )^ / ^ n » {p + qr + {p + q) » n v 
1 {P) 6 J (2n) D ( g + p)2( g 2_ m 2) ^ 9 n-(p + q) 

2 f dD( l 7m (~ 4+ m e ) ( p (p + q) v + (p + q)» n h 



(2ir) D (q + p) 2 (q 2 — ml) \ n ■ (p + q) 

w 2 f d D q 1 



{2ir) D (q + p) 2 (q 2 - m 



2 J (2tt) d (q + p) 2 (q 2 — ml) n ■ q 
We have routed momenta so that — q is the momentum going through the electron prop- 
agator. In the second line of Eq. (|2*H|). we have expanded m = m e — m|/(2m e ). There 
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are three terms in the electron self-energy. The first is the usual electron self-energy in 
n ■ A = gauge. The second term is Lorentz invariant and gives a shift in the electron mass 
of order a m^/m e log (m e //i). Since the integral in the second term is divergent, it gives a 
tiny 0(m 2 p/m 2 e ) correction to the anomalous dimension of the electron mass. The third term 
violates Lorentz invariance, preserves SIM(2) invariance and contains IR divergences due to 
the nonlocality of the theory. We will focus on the evaluation of the last integral in Eq. (|28j). 

Using light-cone coordinates where n- q = q~ and q 2 = q + q~ — q 2 , and restoring the ie's, 
we find 

d D q 1 
(2tt) d n ■ q(q 2 — m 2 + ie)((p + q) 2 + ie) 
d D ~ 2 qdq + dq~ 1 

2(2vr) D q~{q + q- - q 2 - m 2 e + ie)((p~ + q~){p + + q+) - (p + q) 2 + ie) 

where the boldface vectors indicate spatial coordinates orthogonal to both light-cone coor- 
dinates. We evaluate the q + integral using contour integration. The poles in the complex q + 
plane have to lie on opposite sides of the real axis in order for the contour integral in the q + 
plane not to vanish. Therefore, if we close the contour around the pole at q + = (q 2 + m 2 ) / q~, 
the remaining q~ integral is constrained to lie within the range — p~ < q~ < 0, if p~ > 0, 
or < q~ < \p~\, if p~ < 0. Without loss of generality, we assume p~ > and close the 
contour around the pole at q + = (q 2 + m 2 )/q~ . Then we shift q to complete the square in 
the denominator and define the variable, x, 

x = ^-. (29) 

p- 

In terms of this variable, the result can be written as 

*dx r d°- 2 q 1 



n ■ p Jo x J 2(2ir) D 1 q 2 + (1 — x)m 2 — x(l — x)p 2 
i 



l ^(l i _ log ((l-x)ml-x(l-x)p 2 ' ] j (3Q) 



16ir 2 n ■ p J x \e \ fi 2 

where we evaluated the integral over q in dimensional regularization with D = 4 — 2e and 
1/e = 1/e — 7e + log(47r). The x integral is infrared divergent. This divergence is not 
regulated by dimensional regularization and cannot be removed by taking the electron off- 
shell. From our definition oi x, it is clear that the origin of this IR divergence is from the 
limit n ■ q — > 0. In this limit the Lorentz violating term in the electron propagator blows 
up. Rather than using contour integration to evaluate the integral we could have combined 
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the denominators using Feynman parameters, then we would have obtained an expression 
identical to Eq. ([3U|h where x in this evaluation corresponds to Feynman parameter. In this 
way of doing the integral the physical origin of the IR divergence is somewhat obscure. 

The "+" -prescription for defining the SIM (2) invariant theory deals with the IR diver- 
gence by replacing the pole in the electron propagator with the -(—distribution: 



n ■ q 



n ■ q 



(31) 



The integral becomes 
I = 



lQn 2 n ■ p 
i 

16n 2 n ■ p 



dx 
1 dx 



- - log 

e 



(1 — x)m 2 e — x(l — x)p 2 



o 



X 



log I 1 — x — x(l 



X 



p 



mi 



(32) 



so the correction to — zS(p) coming from the third term is 



AE(p) 



a 2 
4tt 13 



n ■ p 



7T 

~6~ 



Li 5 



p- 



mi 



(33) 



Note that now the correction is both UV and IR finite. Going on the mass shell, p 2 
we find the following correction to the SIM(2) invariant electron mass; 

air 



5m 2 p 



6 p 



(34) 



As stated earlier, this correction is finite. 

Next we evaluate the one-loop corrections to the photon propagator. The one-loop dia- 
grams are shown in Fig. ^ the Feynman rules needed to evaluate these diagrams are derived 
in the Appendix. Both diagrams are necessary to obtain a result that satisfies the QED 
Ward identity, p^IP^Qo) = 0. There are 0(mV) corrections to both the vertices and the 
electron propagators. Keeping all terms 0(m 2 g) terms, and denoting the 0(m%) correction 



to the photon self-energy as AIP 1 ' we find 

zAn"" = 



-2e 2 ml 



d D l 



^ n ll p v + p^n v 
® n ■ p 

1 1 



P 



(n ■ p) 2 
n ■ p n ■ p 



x 



(35) 



(2ir) D I 2 — m 2 (I + p) 2 — m 2 e \n ■ I n- (I +p). 

Changing variables / — > — / — p shows that the two integrals in Eq. (j33j) are identical. The 
integrals appearing in the photon polarization are essentially identical to those appearing in 
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FIG. 1: One-loop diagrams for the photon self-energy in VSRQED. 

the electron self-energy and are IR divergent due to the poles inn- 1 and n - (l+p). Note that 
in this case these factors come from both the electron propagator and the SIM(2) invariant 
modifications to the electron photon coupling. Applying the "+" prescription yields: 

n^p" + p^n v p 2 u \ f 1 dx / 



—ml -g» v + — v - ^ y / — log 1 - x{\ - x A 

7r ' \ n-p [n ■ pY J J Q x \ m A t 



This integral can be evaluated in terms of dilogarithms. It clearly vanishes on-shell when 
p 2 = 0. There is no photon mass term generated in C invariant VSRQED at one loop, and 
no Lorentz violating modification of light propagation. 

Divergences of the type encountered in the SIM(2) invariant theory also arise in ordinary 
gauge theory when quantized in n ■ A = gauge. In this case the gauge field propagator has 



poles in l/{n-p) and a prescription is required to define loop integrals 



proposed in Ref. 



22j . The prescription 



is to define loop integrals by making the following substitution: 
1 n ■ q 



n ■ q {n ■ q)(n • q) + ie 

where n M is defined to be (1, — n) if n M = (1, n). This prescription yields sensible results for 
gauge theories quantized in light-cone gauge but will not work in the present case because 
it violates SIM(2) invariance. For example, this prescription leads to the following integral: 

d D q fi i 2p ■ n 1 



(2ir) D (q—p) 2 q-n (Att) 2 n-n e 

The left-hand side is SIM (2) invariant, but the right-hand side is not, because the four- vector 
is not left invariant by all of the generators in the little group of n M . 
The prescription developed here for handling IR divergences is clearly ad-hoc but was 

forced upon us by the necessity of obtaining sensible results from one-loop calculations in the 



18 



SIM(2) invariant theory. Our proposal is only a preliminary attempt to try to address this 
problem. It is satisfactory that one-loop corrections to the photon and electron propagators 
are well-defined, free of IR divergences and for the photon propagator respect the QED Ward 
identity. Interestingly, the corrections to the electron mass is finite, suggesting that if this 
prescription is valid a SIM (2) invariant mass term is technically natural. Several questions 
remain unanswered. Will the "-(-"-prescription continue to yield sensible results in higher 
order diagrams? Will the "-(-"-prescription respect gauge invariance? Though the QED 
Ward identity for the photon propagator is satisfied at one loop, we have not shown that 
Ward identities will be respected in higher order diagrams or for other Green's functions. 
Are there other sensible prescriptions and do they give different physical results? Clearly, 
these issues must be addressed before we can be satisfied that the SIM(2) invariant theory 
is well-defined at the quantum level. 



IV. LEPTON FAMILY NUMBER VIOLATION: fj, e + 7 

In this section we consider the decay /i — > e + 7, which is allowed because the SIM(2) 
invariant mass terms in the lepton sector are not family diagonal. We will not impose C 
invariance but consider the more general Lagrangian in Eq. (|7J). There is a direct fi — > e + 7 
coupling from the covariant derivative in the off-diagonal SIM(2) mass term. This coupling 
vanishes in n-A = gauge. There are also two diagrams which arise due to fi — > e transitions 
induced by the off-diagonal SIM (2) invariant mass term combined with an ordinary QED 
photon coupling to the electron or muon. The Feynman rule for the /1 — > e transition and 
the decay diagrams are shown in Fig. El The factors and in Fig. El are given by 

The left-handed couplings are quite small because they are proportional to differences of 
neutrino masses squared since the off-diagonal mass terms vanish in the limit of degenerate 
neutrinos. For example, using the parametrization of the neutrino mixing matrix in the 
PDG Q, for Mf 2 we find 

M i2 = (VMmlV^)i2 = S12 C12 ci 3 c 23 Am^ + e~* 5 si 3 s 23 c^Am^ - s\ 2 Ara^) 

= -^=si 2 ci2 Am^ + ^=e^si 3 (Am^ - s{ 2 Amy) . (37) 
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FIG. 2: Feynman rule for fermion flavor changing vertex and associated amplitudes for fi~ — * e~j 
at tree level. The third diagram comes from the expansion ^5 term in powers of A and does not 
appear in n ■ A = gauge. 

In the second line we make the approximation S23 ~ C23 ~ l/\^2,Cxs ~ 1- From the cur- 
rent values in the PDG [3], we know that the magnitude of the first term in the last 
line of the Eq. $£ft is 2.6 10" 5 eV 2 while the magnitude of the second term is in the range 
3.0 — 4.710 _4 eV 2 . The relative phase of the two terms is unknown. Assuming that the 
decay is saturated by the contribution from the left-handed SIM(2) mass terms we obtain 
lower bounds for Br[/x — > e + 7] that are 12 to 15 orders of magnitude below the current 
experimental bounds. On the other hand, the right handed couplings are completely unde- 
termined and so the the current bound on Br [/a — > e + 7] gives a nontrivial constraint on 

(UlM 2 U M )i2. 

For the differential decay rate, we find 

dr[/i^e + 7] a f (1 + x 2 )(\M[ 2 \ 2 + |M^| 2 ) - AxRe(M[ 2 M$)\ (n ■ q) 2 



dVL \Q-nm' i pi \ (1 — x 2 ) J n ■ p e n ■ 

where q, p e , and p^ are the momentum of the photon, electron, and muon, respectively. 
We have used the ordinary QED wavefunctions and propagators for the electron and muon 
because the SIM(2) invariant corrections to these give contributions which are suppressed by 
more powers of m 2 , that can be neglected. Lorentz violation manifests itself in the angular 
dependence of the differential decay rate, which vanishes when the photon is emitted along 
the direction defined by n M . Note the result is invariant under rescaling of n M as required 
by SIM(2) invariance. This allows to one to choose n M = (1, n) where n is a unit vector. 
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Performing the angular integral is straightforward and the result for the total decay rate is 



I>-e + 7 ] = (39) 
a (l + x 2 )(|M 1 J 2 | 2 + |Mg| 2 )-4 3 :Re(M 1 j 2 M 1 f) / /J_\ _ 3 2 _ 
Ami (1-x 2 ) 2 V WJ 2 + 2 

where x = m e /m^. Since x ~ 0.005, if \M(^\ < 50 |M^|, then the combination of |M^| 2 — 
4xRe(M ] / 2 M^) will be positive definite. In this case, the contribution from left-handed 
SIM(2) mass terms can be used to establish a lower bound on the branching ratio for 
H — > e + 7, which we find to be 

Br[^ -> e + 7] > 3.2 10~ 26 - 1.0 10" 23 . (40) 

The first number corresponds to S13 = while the second corresponds to S13 equal to its 
current upper bound. These branching ratios are a factor of 10~ 15 -10~ 12 smaller than the 
current experimental bound, Br[/i — > e + 7] < 1.2 10~ n 12[, and therefore uninteresting. 

The smallness of this branching ratio when M[ 2 dominates the decay is directly related 
to the smallness of the neutrino mass squared differences, Amfj. will be proportional 
to the mass squared differences of the eigenvalues of M 2 which are the SIM(2) invariant 
contributions to the right handed lepton masses. Earlier we saw that g e — 2 experiments with 



trapped electrons constrained the flavor diagonal components: (Vm^V^ii ~ (Uj l M 2 UR 
The nonobservation of fi — > e + 7 puts a constraint on the off-diagonal component, M^. We 
find Afg = (U ] R M 2 U R ) l2 < 0.50 10 3 eV 2 . Note that 0.510 3 eV 2 = (22.4eV) 2 , so the upper 
limit on implied by fi — > e + 7 is significantly larger than the bound obtained for the 
flavor diagonal term, M([. 
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V. SUMMARY 



In this paper we have studied the implications of SU(2)l x £7(1) gauge symmetry for 
the SIM(2) invariant neutrino mass proposed in Ref. {oj. This mass term induces Lorentz 
violating effects in QED as well as tree-level lepton family number violating interactions. 
Tests of Lorentz invariance in experiments on trapped electrons severely constrain a possible 
C violating SIM(2) invariant mass term in QED. We studied VSRQED, which contains a 
C invariant SIM (2) invariant mass term for the electron. We derived a Hamiltonian for 
nonrelativistic electrons interacting with background electric and magnetic fields and used 
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it to determine the impact of VSRQED on g e — 2 measurements with trapped electrons as well 
as the spectroscopy of hydrogenic atoms. Corrections from this Hamiltonian are suppressed 
by e = m 2 p/m 2 e . This parameter is expected to lie within the range 9.2 1CT 17 < e < 3.8 1CT 12 
based on our present understanding of neutrino masses. Experiments on trapped electrons 
place limits on e that lie just above the upper limit of this range, while constraints from the 
spectroscopy of Hydrogenic atoms give weaker bounds on e. Bounds on a Lorentz violating 
electric dipole moment comparable in size to existing bounds on the electron electric dipole 
moment could falsify the SIM(2) invariant neutrino mass hypothesis. We found the tree 
level rate for /i — > e + 7 is many orders of magnitude below current experimental limits. 

We also examined one loop self-energy diagrams in VSRQED. These suffer from unregu- 
lated IR divergences due to the nonlocality of the SIM(2) invariant action, indicating that 
the theory is ill-defined at the quantum level. We suggested a prescription for defining loop 
integrals which gives sensible results for the one-loop diagrams considered in this paper. 
QED Ward identities are satisfied, no Lorentz violating effects in the photon propagator are 
induced, and the one-loop correction to the electron's SIM(2) invariant Lorentz violating 
mass is finite. Future work on SIM(2) invariant theories needs to address whether this pre- 
scription, or some other method of defining loop integrals, will give satisfactory results to 
all orders and for all Green's functions. 
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APPENDIX A: FEYNMAN RULES FOR VSRQED 

The Lagrangian for VSRQED is 

where = <9 M — ieA^ is a covariant derivative. There is no modification to the gauge 
field kinetic terms. The parameter nip is (VmitiIV^u = J2i\^ei\ 2m l,i- This is the so- 
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V v 1 — mi + ie 



P 



p z + ie \ n ■ p 



FIG. 3: Light cone gauge Feynman rules in VSRQED. Unlisted rules are unmodified from standard 
QED. 

called electron-based neutrino mass squared. The nonlocal SIM(2) invariant mass term 
gives the correction to the propagator in Eq. (fTrjj) and also gives rise to additional couplings 
of electrons and photons. The additional couplings are absent if the theory is quantized in 
light-cone axial (n ■ A — 0) gauge. In many cases, it is simplest to use light-cone gauge 
when analyzing this theory. Of course, to check that gauge invariance is satisfied we need to 
be able to calculate in arbitrary gauges. We will provide Feynman rules in both light-cone 
gauge and covariant gauges in this Appendix. We also check our covariant gauge Feynman 
rules by verifying that the tree level VSRQED Compton scattering amplitude satisfies the 
Ward identity. 

The propagators for the electron and photon in light-cone gauge are shown in Fig. El and 
the electron-photon coupling is standard. The Feynman rules are simplest in this gauge. 

In covariant gauges, the photon propagator is standard and the electron propagator 
remains the same as in light-cone gauge, but there are now an infinite number of interactions 
due to the term. We can make use of Wilson lines to express these terms in a form that 
is convenient for deriving the Feynman rules. The Wilson line, 

W(x) = exp ^— ie J d\n-A(x + n\) 

satisfies the equation in ■ DW(x) = and is a unitary operator. As an operator equation 
we have 

W\x)n ■ DW(x) = W\x) [n ■ DW(x)} + W\x)W(x)n ■ d 
— n • d , 

which leads to |lf| 

W(x)^-~W\x) . (A2) 



in ■ D in ■ d 
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Expanding the right-hand side in powers of A will give the Feynman rules. In order to 
expand W(x) we need the intermediate integral 

r-0 r rO 



-le 



dX n ■ A{x + nX) = -xe d A pn- A{p) / dX e - ip < x+nX) 

J J -co 

= -ie [ d 4 pn- A(p)e- ip x 



n ■ p + te 

So expanding (jA2|) to first order in A in momentum space (remembering that id in momen- 
tum space will become the sum of the momenta of all operators to its right), we find 

w 1 w \ _^ en-A(q) 1 1 en-A(q) = en-A(q) ( _L 1_ 

in • d n ■ q n ■ p n ■ (p + q) n ■ q n ■ q \n ■ p n ■ p' 

en ■ A(q) n ■ {p' — p) 
n ■ q (n ■ p)(n ■ p') 
e n ■ A(q) 



(n ■ p)(n ■ p') 

Collecting this term with the traditional QED interaction term, the QED vertex is modified 
as shown in Fig. 0] The structure of this vertex is very similar to the modification of the 
weak current derived in the analysis of Ref. In addition to the modification of the 
QED vertex, there are also additional vertices with arbitrary numbers of photon fields. For 
example, second order term in A gives 

e 2 n ■ A( qi )n ■ A(q 2 ) — + 

\n ■ qi n ■ q2 n ■ p \n ■ qi n ■ [p + q2) n ■ q 2 

1 1 1 \ 1 11 



n ■ q 2 n ■ (p + qi) n • qi J n ■ (p + qi + q 2 ) n ■ q x n ■ q 2 

e 2 n ■ A(q 1 )n ■ A(q 2 ) ( 1 ( 1 1^1 



(n-qi)(n-q 2 ) \n ■ p \n ■ (p + qi) n ■ (p + q 2 ) J n ■ p' 

which leads to the two-photon vertex shown in Fig. 0] 

We can check these Feynman rules by verifying that the tree level amplitude for Compton 
scattering obeys the Ward identity. The VSRQED amplitude for Compton scattering at tree 
level in Feynman gauge is a sum of three diagrams shown in Fig. El We define the incoming 
electron and photon momentum to be p and q, respectively, and define the outgoing electron 
and photon momentum to be p' and q', respectively. Next define 



nr 



P H 



2 n ■ p 



(A3) 
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FIG. 4: Covariant gauge Feynman rules in VSR QED up to two powers of A. Unlisted rules are 
unmodified from standard QED. 






FIG. 5: Leading order Compton scattering diagrams in VSRQED. 



and 



2 n ■ kn ■ k' 



(A4) 



Note that 



2 2 2 2 2 

n p — m = p —mp — mQ 



2 2 

p — m e 



(fi p - m )u n (p) = . 



(A5) 
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The first two diagrams in Fig. |3] then contribute 

iM = -ie 2 u n (p) (r u (p + q, p') v 2 +q ^ m \ + q, p) 

V n p+q m 

+r fl (p> - q,p') ^ q + m ° 2 T v {p - q',p)) u n (p) e*»(q')e»(q) . (A6) 

n p'-q ~ m J 

Here & and e* v are photon polarization vectors. The Ward identity demands that the 
amplitude vanish under the replacement e^(q) — *> q** '. Making this substitution in Eq. ()A6|) 
and using the identities 

q^vip + qiP) = ip+q-h 

q^Y^p' - q, p') = - fi p ,_ g , (A7) 



along with the identities in Eq. (jA5|) . it is easy to show that the result is 
iM[e v {q)^q v ] = -ie 2 u n {p')(v v {p + q,p') -T u {jp - q',p)\u n {p) e*" '(<?') (A8) 



2 m l- , ,w ( 1 



-ie —u n (p')jin v — — u n (p) e* u (q f ) . 

2 \n-{p + q)n-p' n ■ {p — q ) n ■ p ' 



The third diagram in Fig. |5] gives the following contribution: 



2 m| _ ^ A fin u ( 1 1 1 1 



iM\e\q) - t?\ = - ie ^u n (p')^ + -— u n (p) e^(q') 

2 n ■ q \n ■ p n ■ p' n- (p + q ) n - p — q ' 



ie 2 ^u n {p') fin u u n (p) e w (q') , (A9) 

2 \n-{p + q)n-p' n ■ [p — q')n ■ p ' 



which exactly cancels Eq. (|A8|) . so the Ward identity is satisfied. 
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